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ABSTRACT 

Thanks to a stellar evolution code able to compute through the C-flash we link the 
binary population synthesis of single degenerate progenitors of type la supernovae 
(SNe la) to their physical condition at the time of ignition. We show that there is a 
large range of possible ignition densities and we detail how their probability distribu- 
tion depends on the accretion properties. The low density peak of this distribution 
qualitatively reminds of the clustering of the luminosities of Branch-normal SNe la. 
We tighten the possible range of initial physical conditions for explosion models: they 
form a one-parameter family, independent of the metallicity. We discuss how these 
results may be modified if we were to relax our hypothesis of a permanent Hachisu 
wind or if we were to include electron captures. 
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1 INTRODUCTION 

Phillips relations iPhil hpslll993l) have made type la super- 
novae (SNe la) a very useful tool for cosmology. The corre- 
lation between their absolute peak brightness {M pca ^) and 
their decline rate (Amis) has been used to infer distances of 
high redshift galaxies and refine our know ledge of cosmolog- 
ical parameters (e.g. iTonrv et al. lEoOlT) . However current 
models are not able to reproduce the range of Ami5 ob- 
served: we still do not understand what primary parameter 
is responsible for the main diversity of SNe la. Moreover 
there exists a finite spread of M pca k for a given Amis which 
hints at the existence of a secondary parameter. Knowledge 
of both these controlling parameters could help remove any 
bias from cosmological measurements due to an evolution of 
the secondary parameter and help reduce the scatter of the 
Phillips relation. 

With synthetic light curves and spectra with explosion 
models, attempts have been made to infer what parame- 
ters of the explosion contro l which observational property 
(Rookc & Hillcbrandt 2004, for example). In this paper we 
proceed from the binary population synthesis of SNe la pro- 
genitors where we only consider the single degenerate chan- 
nel, one of the most promising channels for the progenitors of 
SNe la. We model their evolution until we get the physical 
conditions at the time of the ignition. We start our com- 
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putations just after a white dwarf (WD) has been formed 
through a common envelope evolution phase. The initial pa- 
rameters are the WD composition, its initial mass Af WD , the 
initial mass M2 of the secondary and the initial separation 
cii. The WD cools down while the secondary evolves until 
it overflows its Roche lobe at the end of its main-sequence 
evolution (cooling phase). Accretion on to the WD heats it 
up and enhances its central density until the carbon fusion 
starts (beginning of the accretion phase) . The rate of energy 
production due to burning and accretion soon overcomes the 
losses due to escaping neutrinos and photons: thermal bal- 
ance does not hold anymore (birth of the convective core). 
However, the high electron degeneracy and the high temper- 
ature sensitivity of the burning rate prevent expansion from 
controlling the burning and the convective core has to grow 
very fast to cope with the increasing energy generation rate: 
a flash ensues (C-flash phase). Convective flows are able to 
get rid of the energy released by carbon burning until the 
turnover time scales become too long relative to the heating 
rate. At this point it is believed that one or multiple bub- 
bles ignite and a flame front propagates out and eventually 
unbinds the star (ignition or beginning of the explosion). In 
the present work we investigate the probability distributions 
and correlations between various characteristics of the star 
at the time of ignition which are thought to influence the 
outcome of the explosion. 

Section 2 presents the details of our numerical setup and 
the main grid of models. In section 3 we describe the result- 
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ing ignition conditions. In sections 4 and 5 we investigate the 
role of metallicity and the effects of binary evolution respec- 
tively. And in sections 6 and 7 we discuss and summarise 
our results. 



2 METHOD 

2.1 Stellar evolution code 

We im proved the Eggleton code jEggletorJIlimllPols et al~l 

I1995D to enable the modelling of white dwarfs from their 
cooling tracks and the start of the accretion to the very late 
phases of the C-flash until shortly before the explosion takes 
place . Note that neutr ino losses in the Eggleton code are 
from lltoh et al. I lll992l) : they are an important factor that 
determine the position of the central density-temperature 
tracks of the WD. 

We consider only the 12 C+ 12 C fusion reaction and as- 
sume that 24 Mg is its only product. We hence completely ne- 
glect all nucleosynthesis issues such as Urea processes which 
are l ikely to occur during the C-flash (see ILesaffre et al. I 
120051) . 

For each chemical species described in the code we solve 
not only for its mass fraction but also for its spatial gradient. 
This allows us to compute the mixing for realistic diffusion 
time scales as given by the mixing length theory (MLT). This 
is important because close to ignition convective turnover 
time scales become comparable and even longer than the 
burning time scales. 

We compute the equations of hydrostatic stellar evo- 
lution on a staggered mesh which allows us to take arbi- 
trarily small time steps required by the burning time scales. 
This is also needed to stabilise chemical equations including 
both gradients and composition. Radius, mass, luminosity 
and chemical gradients are hence computed on the inter- 
faces of the computational shells whereas all other variables 
are taken at the centre of mass of each zone. 

We make use of a movin g grid algorithm JPorfi fc Drurvl 
Il987l: ILesaffre et al. lEo04ah with a time delay of one year. 
This friction of the grid limits numerical diffusion and ap- 
pears to stabilise the code against rapid changes at the onset 
of accretion and near ignition. 

Last but not least we implemented a scheme to track 
the boundaries of convective regions within a zone. As the 
time scale for burning decreases, a discrepancy builds up be- 
tween the rate of change of temperature within and outside 
the convective core. When convective boundaries are tied up 
to the interfaces of the computational zones, the code has 
more and more difficulty to converge and breaks down way 
before the explosion takes place. We need to release the con- 
vective boundaries and interpolate the convective criterion 1 
to infer their position inside a zone. We are then able to ac- 
count properly for the discontinuity in the rate of change of 
temperature. Chemical rates of change benefit from exactly 
the same treatment but chemical profiles are discontinuous 
at convective frontiers (unlike the temperature) and they re- 
quire an additional flux through the advancing boundary of 
the core. 



We use the Schwarzschild criterion without overshooting 
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Figure 1. Age vs. mass of the secondary for our sample of type 
la progenitors. The solid and dashed lines are analytical esti- 
mates for the age of the main sequence t aken from lEggleton et al. I 
il98fll) (solid) and lHurlev et al. I i200(l (dashed). 

2.2 Initial WD models 

We build a series of WD models of different masses M^yo by 
accreting on a very low-mass (0.3 M0) C + O WD model at a 
rate of 10 -7 Mq /yr which mimics the growth of the CO core 
during the giant branch evolution. This gives a very crude 
estimate for the thermal state of the WD when it is first 
formed. This is however not essential as the thermal profile 
relaxes duri ng the subsequent cooling phase. lUmeda et al. I 
(1999a) and iDomfnguez et al. I i200lT) have shown how the 
mass and metallicity of the progenitor star of the WD in- 
fluence its chemical state. In the present study we use their 
results to build synthetic chemical profiles of the initial white 
dwarf. We use their central values for the C/O ratio (linearly 
interpolated at My^) to build a central core of mass half 
the white dwarf. This is then connected linearly to a C/O 
ratio of 1 uniform in the outermost 10% in mass of the WD. 

2.3 Cooling phase 

The newly formed WD cools down at constant mass until the 
secondary star fills its Roche lobe, which generally happens 
at the end of its main sequence, on its way to become a 
giant (in the Hertzsprung gap) or on the giant branch. In 
each case the cooling time is close to the main-sequence 
age of the secondary star (see figure . For a given WD 
mass the length of the cooling phase generally determines 
the temperature at its centre. It does not depend much on 
its initial thermal profile. Hence the initial conditions for the 
accretion phase depend only on the age of the secondary and 
the mass of the WD. 

2.4 Accretion phase 

When the secondary starts to overflow its Roche lobe 
it transfers mass on to the WD. How fast an accretion 
rate the WD can cope with is a crucial question and is 
still the object of m u ch de bate. In this study we follow 
lHan fc Podsiadlowskil §004) and assume that the Hachisu 
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wind model jHachisu et al. i ll 996) holds. The WD accom- 
modates high mass transfer rates by getting rid of the extra 
mass through a wind as long as the mass transfer rate ex- 
ceeds a critical rate M CT given by 



5.3 x 10" 



■ (1-7-X) 
X 



(M WD /M - 0.4) M yr" 1 (1) 



where X — 0.7 is the hydrogen mass fraction in the accreted 
material from the secondary and A/wd is the mass of the 
WD. As a further simplification we assume that the mass 
transfer rate is always higher than this critical rate so that 
the Hachisu wind determines the accretion rate on to the 
WD. The net growth of the C+O WD is finally modulated 
by the efficiency of He-shell flashes according to 

MwD = ??He(Mor) x M cr (2) 

where rjHe is given by equation (4) of lHan fc Podsiadlowski 
i2004h . With these assumptions the growth rate of the WD 
depends only on the mass of the WD. We estimate the effects 
of the end of the wind phase in section 5. 



2.5 Ignition 

When convection is no longer able to evacuate the energy 
input from the C-burning, thermal runaway starts in some 
bubbles near the centre and the star explodes. Exactly how 
and when this occurs is still unclear but recent developments 
allow us t o pinpoint analytic conditions for the ignition o f 
the flame JWooslev et al. 112004: IWunsch fc Wooslevll2004D . 
Here, we stop our computation when the differential burn- 
ing time scale t\, is shorter by a given fraction a of the local 
convective turnover time scale t c at some point in the star. 
Here th is the e-folding time scale of a temperature fluctua- 
tion subject to C-burning: 



1 



dq 



th ~ cpTdlnT^p' (3) 

where T is the temperature, p is the mass density, cp is the 
specific heat at constant pressure and q is the rate of energy 
generation due to C-burning and t c is a convective element's 
crossing time over a pressure scale height: 

Hp 



tc = 



(4) 



where Hp is the pressure scale height and u c is the convective 
velocity as given by MLT. 

When th — t c (a — 1) the differential burning time scale 
is of the same order as the convective turnover time scale. 
At this point differential reactivity (the fact that upward 
and downward moving fluid elements burn at different rates) 
should be included in the treatment of conve ction to get ac- 
curate mixin g and temperature excesses (see lLesaffre et al~l 
l2005l.l2004bh . 

In our models (in which we use classical MLT) it turns 
out that the temperature excess is typically 0.025% when 
this happens. It takes roughly eight e-foldings to grow from 
0.025% to 100%. Hence can we estimate that temperature 
fluctuations will be big for % = t c /8 (a = 1/8). At this 
point, a typical ascending temperature fluctuation increases 
its temperature up to twice the average temperature before 
it is mixed. We take this as a plausible point for the start of 
the thermal runaway. 



IWunsch fc Woosievl i2004l) consider a temperature 
background fixed in time and analyse the growth of per- 
turbations relative to this background. As a result they find 
that the relevant time scale for the growth of the temper- 
ature fluctuations is the e-folding time scale for the overall 
burning rate 



t- 1 - q 



c P T 



(5) 



Because of the high power exponent in the temperature de- 
pendence of the burning rate th — trj/22. Hence in this 
picture a = 1/22 is a more relevant requirement for the 
start of the explosion. Note that tp, also indicates the time 
scale of variation for the average temperature in the cen- 
tre of the star. When tp, < t c a time-dependent convective 
model should be used because the turnover time scales be- 
come longer than the global evolutionary time scales. 

To account for the uncertainty in what determines the 
start of the explosion, we record models when a=l,l/8 and 
1/22 as possible times for the start of the explosion. In prac- 
tice all our simulations reach the a = 1 point, half of them 
reach a = 1/8 but all break down before a — 1/22. We 
hence extrapolate our results from our last converged model 
to get the conditions for all these different possibilities for 
the ignition point. In almost all our simulations the criterion 
for explosion is first fulfilled at the centre because the local 
convective time scale is longer there. The only exception is 
for one run which started off-centre carbon burning and the 
convective zone never reached the centre (its parameters are 
t a = 0.8 Gyr and M^d = 0.9 Mq, the bottom of the con- 
vective region is at 0.05 Mq at ignition). We display results 
for this run only in figure 0] as most of the quantities that 
take part in the correlations are defined at or relative to the 
centre. 



2.6 Parameter space 

Because we always assume a Hachisu wind phase (except in 
section 5), the entire evolution of the WD is determined by 
the cooling age and the initial mass of the WD. Starting with 
a given initial model, we let it cool down for a time t a after 
which we trigger the accretion rate Mwd- The beginning of 
the accretion is smoothed linearly w ith a short time scale 
(10 5 y r) to mimic the rise observed in lHan fc Podsiadlowski 
(2004). We build a grid of 24 models spanning six differ- 
ent masses M^d evenly distributed in the range [0.7,1.2] x 
M and i a =0.1, 0.2, 0.4 and 0.8 Gyr. This grid of models is 
designed to s ample the parameter space o f SNe la progeni- 
tors found by lHan fc Podsiadlowskil j2004l) . The metallicity 
is Z = 0.02 and we use initial C/O ratios from lUmeda et ahl 
(1999a) for this metallicity. 



3 IGNITION CONDITIONS 

Parameters which have been found to influence the explo- 
sion outcome are density, temperature, carbon mass fraction 
and convective properties. We investigate here how they are 
related to the parameters M^d an d t a - 



4 P. Lesaffre, Z. Han, C. A. Tout, P. Podsiadlowski & R. G. Martin 




Figure 2. Central density-temperature tracks for different initial 
WD masses (from left to right: 0.7, 0.8, 0.9, 1.0, 1.1 and 1.2 Mq) 
and the same cooling age t a = 0.4 Gyr. Squares, triangles, crosses, 
stars and diamonds indicate respectively the start of the accretion 
phase, the growth of the convective core and the a = 1, 1/8 and 
1/22 ignition points. 



Figure 4. Ignition densities (denned by a = 1) vs. initial 
WD mass for different cooling ages t a = 10 s yr (diamonds), 
t a = 0.2 Gyr (triangles), t a = 0.4 Gyr (squares) and t a = 0.8 Gyr 
(crosses). Densities are only slightly lower for other values of a. 
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T c track because the accretion rate depends only on A/wd- 
As a result, they all have the same ignition density. This is 
reflected in the distribution of ignition densities which has 
a spike at the lowest density. 

Densities also seem to all converge towards the same 
maximum density for higher initial masses M^yrj. This is 
due to the shape of the ignition curve which is almost verti- 
cal in the p — T plane because electron screening enhances 
C-burning at high density. As a result, the p c — T c tracks that 
go to high densities tend to all start the C-flash phase at the 
same density. Note also that these high densities point to- 
wards the importance of electron captures, not yet included 
in this study. 

Ignition density is the parameter that varies most along 
our grid of parameters. We shall show that most of the other 
physical properties at ignition are well correlated with this 
parameter. 



Figure 3. Central density-temperature tracks for different cool- 
ing ages (from left to right: 0.1, 0.2, 0.4 and 0.8 Gyr) and the 
same initial WD mass 1 Mq. Symbols are as in figure 121 

3.1 Central density 

Larger initial masses M^v D yield higher ignition densities be- 
cause their initial densities are higher (see figure|2l . A longer 
t a yields a lower initial temperature and a hence higher ig- 
nition density (see figure 01. 

The dependence of density on Mwd and f a is summa- 
rized in figure 2] A WD with a smaller Mwd has a shorter 
thermal diffusion time scale and a longer accretion phase. 
As a result heat waves have time to propagate through the 
star and the thermal profile of the star reflects the global 
balance between the energy input from the accretion (and 
burning) and the energy losses from neutrinos. In the rest of 
the paper we refer to this state as global thermal balance. 
WDs in a state of global balance all end up on the same p c - 



3.2 Ignition Temperature 

A lower a yields hotter and only slightly less dense ignition 
conditions because matter is still degenerate at this point. 
We now examine how the temperature varies according to 
density when we define the ignition conditions following a 
given a. 

The convective turnover time scale t c is rather indepen- 
dent of any other parameter at the time of ignition as shown 
in section HOI By contrast the electron screening corrections 
to the rate of C burning increase it at high density. As a 
result, because the ignition temperature is defined by the 
relation ib = at c , ignition temperatures are lower for higher 
densities (see figure [KJ. This m ight introduce in t eresti ng ef- 
fects in light of the results of llapichino et al. I (I2005TI who 
find a strong sensitivity to the temperature for the proper- 
ties of igniting bubbles. 
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Figure 5. Central temperature vs. central density at ignition 
central density for a = 1 (solid), 1/8 (dotted) and 1/22 (dashed). 
Curves and signs are otherwise labelled as in figure ^ 
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Figure 6. Central convective time scale t c vs. central density at 
ignition, labels as in figure IB1 



3.3 Convection 

Explosion models have recently focused on the im- 
portance of t he convective state at the beginning of 
the explosion jReinecke et al. 1120021: iHoflich fc Steinll2002l; 
iGarci'a-Senz fe Bravoll2005l : IWunsch fc Wooslevll2004l) . It is 
interesting to note that our simulations (based on MLT) 
show very little variation for the convective state at the time 
of ignition. The turnover time scale t c appears to be more 
or less independent of t a and Afw D . Naturally these time 
scales depend upon our definition for the ignition as t c di- 
rectly enters our criterion (see figure Ej. We recover values in 
the range 10-100 s as quoted bv I Wunsch fc Wooslevl d2004T) 
when we use a = 1/22. Also worth noting is a weak depen- 
dence of t c on central density. This could be enhanced by 
convective Urea processes due to th e feedback of electro n 
captures on convective velocities (see lLesaffre et al~1 [2005 1 . 
The convective core mass at ignition is mainly deter- 
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Figure 7. Convective core mass at ignition vs. central density at 
ignition, labels as in figure l5l 

mined by the central density which sets the pressure-mass 
profile and the central temperature which determines the 
extent of the convective core in the temperature-pressure 
profile. It is hence not surprising that we recover the same 
kind of correlations as for the central temperature (see fig- 
urc0 . We investigate this para meter in light of the rot ating 
progenitor models computed bv lYoon fc Lan gcr (2004) who 
found super- Chandrasekhar solutions with a strongly differ- 
entially rotating layer around the mass shells m = 1.1 Mq 
to m = 1.25 Mq. If the convective core is able to reach this 
shell before ignition, it might suppress the centrifugal sup- 
port and leave a higher than Chandrasekhar mass object 
which may then collapse instead of exploding. On the other 
hand convective instability might be suppressed by differen- 
tial rotation and that would trigger the thermonuclear run- 
away sooner. Figure[7|suggests that the convective core may 
sometimes reach the differentially rotating layer, depending 
on what definition we use for the ignition point. 

3.4 Carbon mass fraction 

The evolution of the central carbon mass fraction is mainly 
driven by the growth of the convective core which mixes 
material richer in carbon at a higher rate than the car- 
bon is consumed at the centre. As a result the final carbon 
mass fraction depends essentially on the composition pro- 
file and the initial mass of the WD, that is on the initial 
Mwd ( see figure The fact that there is little variation 
among different WDs and the recent result that the explo- 
sion might be less sensitive t han we think to this parameter 
jRopke fc HillebranddlitJOl) lead us to conclude that this 
might not be an essential parameter for understanding the 
observational properties of SNe la. However it is the only 
parameter that slightly decouples from the central density. 

3.5 Magnesium mass fraction 

Thanks to our oversimplified treatment of C-burning, the 
24 Mg mass fraction represents the amount of carbon burnt 
during the C-flash. We show that this is weakly dependent 
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Figure 8. Carbon central mass fraction at ignition vs. initial WD 
mass, labels as in figure El 



Figure 10. Histogram of the central density at ignition for our 
5576 SNe la progenitors. 
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Figure 9. Correlation between Magnesium mass fraction and 
central density at ignition, labels as in figure 151 



on the central density even though the trend is not clear at 
high densities (see figure |5J. Later ignition (lower a) natu- 
rally burns more carbon. The amount of carbon burnt is of 
interest with regard to Urea species which are among the 
products of C-burning ( 23 Na and 25 Mg, for example). 



3.6 Distributions 

We use the results of a binary population synthesis that pro- 
duced 5576 SNe la progenitors out of 10 7 initial s ystems with 
«ce = a t h = 1 (see lHan fc Podsiadlowskill2004l section 4). 
Note that the WD+red giant branch channel is not included 
in their study: this may affect the distributions. We take the 
age when the WD reaches the Chandrasekhar mass and the 
initial mass of the WD as parameters (t a ,M^ VD ) which we 
interpolate linearly on our grid of models to get the central 
ignition density pi for each of these systems. We display the 
histogram of pi in figure [TU1 

The very sharp peak at the lowest density reflects the 
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Figure 11. Histogram of the central density at ignition for ages 
greater than 0.8 Gyr (solid), ages between 0.4 and 0.8 Gyr (dot- 
ted) and ages lower than 0.4 Gyr (dashed). 



fact that accretion proceeds on longer than thermal diffusion 
time scales for the lowest initial masses M^n- The secondary 
peak at high density reflects the shape of the ignition curve 
due to electron screening at high density. It is tempting to 
associate each one of these peaks to the clustering of either 
Branch-normal SNe la or 91bg like events. But computations 
of the outcome of the explosions for each of our models need 
to be performed before we can assess the relation between 
ignition densities and observational properties of the explo- 
sion. The shape of the distribution of MJv D is rather flat 
and does not influence the o verall probability distr i butio n 
of pi. However, as noticed by Han & Podsiadlowski (2004), 
is slightly correlated with f a with lower ages giving 
rise to higher WD masses on average. As a result more rel- 
ative weight is given to the high density peak for younger 
systems (see figure [TTt . 
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Figure 12. C/O ratio for lUmeda et al. Ul999al). Z = 0. 02 (solid) 
and Z=0.001 (dotted) and for IppiTiinguez et al. I feOOfl . Z=0.02 
(dashed). Diamonds are the values computed by these authors, 
'+' signs are our interpolated values to compute the initial mod- 
els. 



Figure 13. Central carbon mass fraction at ignition for each 
initial model of figure 1121 Diamonds (solid line) are our refer- 
ence run (see fi gure IS } . triangles (dotted line) are for Z=0.001 
Umcd a et al. 1999a) and squares (dashed line) are Z=0.02 
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4 COMPOSITION EFFECTS 

In this section we discuss the effects of the composition of 
the white dwarfs. The initial mass and initial metallicity of 
their progenitors can vary independently but there can be a 
correlation between initial mass of the progenitor and that 
of the white dwarf at the end of the cooling phase. 



4.1 Initial C/O ratio 

In the previous section all our results are based on cen- 
tral C/O ratios taken from lUrneda et al. I lll999al) . How- 
ever iDommguc z et al. ] fcOOll) have found significantly dif- 
ferent values for the same solar metallicity as shown in 
figure ^| This is probably due t o their use o f differ- 
ent 12 C(q,7) 16 reaction rate s withlUmeda et al~l (Il999a|) 
using 1.5 times the rate of JCaughl an fc Fowler I ll988l) 
whereasjDmnmgiiezetji ] feOOll) use a higher rate given 
bv lCaughl£ni^t^0^n^85ir We co mpute here a grid of mod - 
els using central C/O values from iDommguez' et al. I (I2001D 
for a fixed cooling age of t a = 0.4 Gyr and the same range 
of initial masses as previously. 

The resulting physical conditions at ignition are virtu- 
ally identical to the runs described in the previous section. 
The only slight change is for the central carbon mass frac- 
tion which changes as shown in figure [T"3*l This change only 
reflects the mixing in the convective core of the initial com- 
position profile with the accreted matter. 



4.2 Metallicity effects 

As shown in figures 1121 and 1131 a change in metallicity is 
very much like a change in the initial C/O ratio, except 
that the magnitude of the effect is even less. In this respect, 



evolution in the C-flash and computed as in 


Timmes et al. 1 


( 2003!). However, as noted bvlLesaffre et al. 1 


2005), electron 



C-flash. Whether or not the initial amount of 22 Ne has an 
effect on the detailed nucleosynthesis during the C-flash can 
be addressed with a more extensive nuclear network than 
used in this study. 



5 BINARY EVOLUTION EFFECTS 

The results of previous sections strongly rely on the assump- 
tion of a permanent Hachisu wind, which allows us to neglect 
the binary evolution during the accretion process. Here we 
examine the effect of varying the accretion rate history. 



5.1 End of the wind phase 

From figure 1 of lHan fc Podsiadlowskil (I2004T) we find that 
the mass transfer rate decays exponentially after the wind 
has stopped. We hence model the mass transfer rate as 

M W d - Mo 



M = M cr (M ) 



to 



(6) 



captures could change the neutronisation in the course of the 



where Mo is the mass of the WD when the wind stops and 
to is the time scale for the decay. The value of to varies 
fr om 2.7 x 10 5 yr to 9 x 1 s yr in different panels of figure 1 
of lHan fc Podsiadlowskil (120041) and we take 6 x 10 5 yr as a 
characteristic value. We then use A/wd = r^HcM. 

For three different masses M^d = 0.8, 1.0 and 1.2 Mq 
and t a =0.8 Gyr we compute the evolution after the end of 
the wind phase for Mo = 1.1, 1.2 and 1.3 Mq (except for 
Mwo = 1-2 M for which M = 1-25, 1.3 and 1.35 Mq). 
Only the central density and the convective core mass at 
ignition are found to change significantly compared to the 
models of section 2. 

Figure tTH sums up the results for the p c — T c tracks. The 
slope of a given track in the p c — T c diagram depends on the 
ratio of the rate of change of density and temperature. For 
a degenerate WD the density is directly linked to the total 
mass and hence the central density reacts immediately to 
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Figure 14. Central density-temperature evolution for three dif- 
ferent initial WD masses: 0.8 (left), 1.0 (middle) and 1.2 Mq 
(right). Solid lines are the reference runs with a Hachisu wind al- 
ways on. Other lines stop the wind when the mass reaches three 
different critical values of Mo (see text). Labels are dotted, dashed 
and dash-dotted for decreasing values of Mo ■ 

any change in mass. On the other hand, the central temper- 
ature reacts to a change in accretion at the edge of the WD 
with a delay of the order of the thermal diffusion time scale 
td across the WD. When the accretion rate slows down, the 
rate of change of the central density slows down immediately 
while the temperature continues to increase at the same rate: 
this steepens the slope of the p c — T c track. For the lowest 
masses (left in figure 1141 , the accretion takes place over a 
longer time than td- Hence global thermal balance is soon 
reached again and because the accretion rate is smaller this 
leads to a lower temperature for a given central density. The 
final result is a slightly higher ignition density. For interme- 
diate masses (middle in figure [TH . global thermal balance 
was not realised at all because the accretion rate was too 
fast. Hence the temperature was below the temperature the 
star would reach at global thermal balance. Slowing down 
the rate now gives time for the star to adjust and its tem- 
perature increases, which leads to lower central densities at 
ignition. For the highest masses (right in figure [Lit , there is 
not enough time for the star to react: the effects of slowing 
down the accretion rate do not reach the central thermal 
state of the star. 

A more significant difference exists for the mass of the 
convective core at the time of ignition. Indeed, the lower ac- 
cretion rates yield lower temperatures at the edge of the core. 
Because the outer boundary of the convective core is deter- 
mined by the point where the inner adiabatic temperature 
profile meets the outer temperature profile, the convective 
core is bigger for lower accretion rates at a given central tem- 
perature. We find an increase of between 0.05 and 0.2 Mq 
in the convective core mass at ignition when the end of the 
wind phase is modelled. 

We also computed the evolution of a WD using Mo = 
1.12 Mq and to = 9.15 x 10 5 yr which gives a very accu- 
rate fit to the mass transfer rate for panel (a) of figure 1 of 
lHan fc Podsiadlowskil i2004h . But modelling the end of the 



wind phase actually makes only very little difference in this 
case due to the low initial mass of the WD (0.75 Mq). The 
only significant change was an increase of about 0.07 M© in 
the convective core mass at ignition. 

5.2 Varying the critical accretion rate 

To account for uncertainties in the critical accretion 
rate in the Hachisu wind model (see the discussion in 
lHan fc PodsiadlowskfeOOflF) . we ran the sequence of masses 
for t a =0.4 Gyr with twice the value provided by equation 
The effect turns out to be quite small, with less than a 
10% variation. WDs of lowest initial masses have a slightly 
lower pi because the global thermal balance for a higher 
energy input rate yields higher temperatures. The effect is 
reversed at high initial masses because the shorter accretion 
times leave less time for the star to reach global thermal 
balance, so that the rise in central temperature is even more 
delayed and higher densities are reached at ignition. 

5.3 Qualitative results 

To properly include the end of the wind would require mod- 
elling the stellar evolution of the secondary along with the 
WD. We do not do this explicitly here. However, a decrease 
of the accretion rate at the end of the wind phase as in 
section 15.11 mainly lowers the ignition density for WDs of 
intermediate initial masses. This suggests that binary evo- 
lution effects depopulate intermediate ignition densities to 
the benefit of the lowest densities. Hence the probability dis- 
tribution function of ignition densities should have a sharper 
peak at high density and a broader peak at low density. The 
broadening of the low-density peak would also be empha- 
sised by the slight increase in pi for lowest initial WD masses 
in systems that terminate their wind early on. 



6 DISCUSSION 

The main uncertainty in our model comes from the Hachisu 
wind model (discussed in section |SJ but our results seem to 
hold even for significant variation in the accretion history. 
In section 4 we showed that, although the initial C/O ratios 
are not well known, they do not affect our results very much. 

The main caveat of our study is that we do not include 
electron captures while our models go to quite high densities. 
In the future we hope to be able to include the convective 
feedback as in iLesaffre et al. I (2005) but until now tech ni- 
cal difficulties have prevented us iLesaffre et al~ll2004bl) . If 
there was indeed a negative feedback on convective veloc- 
ities, the mixing would be less efficient. Consequently the 
abundance changes would be enhanced and the explosion 
would probably take place sooner, at higher densities. In par- 
ticular electron captures would increase the neutronisation 
during the C-flash which may h ave dramatic consequ ences 
for the outcome of the explosion llTimmes et al. l2003tl . Less 
homogeneous abundance profiles at the time of the explosion 
may also have consequences for the way the flame propagates 
because it goes through a medium increasingly concentrated 
in fresh fuel. Finally electron captures would tend to slightly 
increase the density compared to our models but this would 
probably be a small effect. 
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The high densities encountered in our models raise 
the question of whether the star will undergo an electron- 
capture supernova type of event leading to the formation of 
a ne utron star rather t han a thermonuclear explosion. How- 
ever [&itferee^^d] j2005t) indicate that a thermonuclear 
explosion is very likely given the high C mass fractions in 
our models. 

Accreting WDs a re almost certainly strongly rotating. 
lYoon fc LangeJ i2004f) have shown how much influence this 
could have on the fate of the WD and how this could account 
for part of the type la diversity. However rotation is directly 
linked to the accretion so that it probably does not generate 
a new independent parameter. Our study has shown, at least 
in some cases, that the convective core is likely to reach the 
strongly differentially rotating region that helps support the 
star to higher than Chandrasekhar masses. Rotation might 
then help to increase the influence of binary evolution effects 
that seem secondary in this study. 

Finally this study is valid for the single degenerate chan- 
nel only. If a Hachisu wind holds as well for the accretion of 
a WD on to another WD, we may apply part of the results 
presented here but in this case the C/O ratio of the accreted 
matter can be different from 1. However the complete dis- 
ruption of the accreted WD is likely to lead to a picture 
quite different from the Hachisu wind. 



7 CONCLUSIONS 

We link the properties of single degenerate type la progeni- 
tors at the time of the ignition to results of a binary popu- 
lation synthesis thanks to successful models of the C-flash. 

We show that there is a large range of possible igni- 
tion densities whose distribution reflects the properties of 
the accretion on to the WD. Furthermore, almost any other 
property of the WD at the time of the ignition is correlated 
to the ignition density: this makes it a one parameter family 
of models which considerably tightens the range of possible 
initial conditions for explosion models. The convective time 
scales, amount of C-burning and central temperatures are 
functions of the central density at ignition. The central C 
mass fraction is better correlated to the initial mass of the 
WD progenitor. The core mass at ignition can be modified 
according to how the mass transfer proceeds (this implies a 
sensitivity to the initial separation of the binary) . These last 
two parameters may hence be regarded as slightly indepen- 
dent of the ignition densities, although they do not have a 
large range of variation and are hence probably less essential 
than the density for the explosion. 

We have investigated several criteria that determine 
when the star would explode. Although the actual depen- 
dence on density of the temperature, convective state and 
amount of C-burning are sensitive to the criterion chosen, 
we note that their correlations to density are still tight and 
the distribution of ignition densities does not change. 

We show that the metallicity has almost no effect on the 
C-flash phase. Hence this parameter is independent of the 
central density at ignition. Because both have been shown 
to have a big effect on the explosion outcome, we postu- 
late that ignition density and metallicity are the two main 
parameters responsible for the diversity of type la super- 
novae. Metallicity has a linear effect on the peak bright- 



ness of type l a supernovae with low metallicity yielding 
brighter M„ c „i, jTimmes et al. lEool iTravaglio et al. 112003 ; 
iMazzali fc Podsiadlowski l2005h . However the effect of den- 
sity on Mp 0a k is still unclear: hydrodynamics of the explosion 
suggest that more mass is burnt to nuclear statistical equi- 
librium for higher den sities but explosion nucleosynthesis 
jTravaglio et al. 112004) predicts that more stable elements 
are produced at higher densities because of electron cap- 
tures. 

If density is the primary parameter and metallicity is 
the secondary parameter, then metallicity changes could 
shift the Phillips relation. On the other hand if metallicity 
is the primary parameter and density is the secondary pa- 
rameter there should be less change in the Phillips relations, 
unless metalli city strongly af f ects the M^d distribution as 
postulated bv lUmeda et al.1 lll999bl) . Finally, density and 
metallicity might be degenerate and both be primary pa- 
rameters. The chain of computations from the ignition con- 
ditions to the light curves via explosion models and explosive 
nucleosynthesis (including electron captures) must be com- 
pleted before we can theoretically address the effect of each 
parameter and uncover the physical nature of the Phillips re- 
lation. Observations of correlat ions with host galaxy metal- 
licities JCallagher et al. 1120051) and comparison of type la 
luminosity functions with the metallicity or ignition density 
distributions will also give important clues. 
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